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Abstract 

A modified gravitational theory is developed in which the gravitational coupling constants G and Q 
and the effective mass m,/, of a repulsive vector field run with momentum scale k or length scale l = 1/As, 
according to a renormalization group flow. The theory can explain cosmological early universe data with 
a dark hidden photon and late time galaxy and cluster dynamics without dark matter. The theory agrees 
with solar system and binary pulsar observations. 

PACS: 04.60.Be; 04.60.Gw 


1 Introduction 

A modified gravity theory called Scalar-Tensor-Vector-Gravity STVG II] has been developed in which the 
gravitational coupling strength G, and the mass-range g of a vector field q were promoted to dynamical 
scalar fields in the action. The gravitational charge of the repulsive vector field is Q = \ZchGnM, where a 
is a parameter, Gn is Newton’s gravitational constant and M is the mass. The structure growth and CMB 
data in early universe cosmology were explained by postulating a hidden, dark photon 7 ' identified with 
the spin 1 vector field that plays the role of cold dark matter, non-thermally produced in the very early 
universe [ 2 ]. As the epoch of the formation of the first galaxies evolved the effective mass y decreased in 
size and its late time value was fixed by fits to galaxy rotation curves and cluster dynamics without dark 
matter urns to be n = 0.042 kpc , corresponding to a mass = 2.6 x 10 28 eV. The bullet cluster 

dynamics has been explained without dark matter | 6 j. By promoting G and y to be dynamical scalar fields 
in the action and field equations, it was necessary to postulate screening mechanisms to guarantee agreement 
with solar system and binary pulsar timing data [7|. 

In the standard model of particle physics, the coupling constants cm = g 2 / 47 r run with momentum scale k 
according to a renormalization group (RG) flow scenario based on a Wilsonian action. The coupling constants 
ai{k) are not promoted to be dynamical fields in the Abelian and Yang-Mills actions. For the non-Abelian 
fields the screening leads to asymptotic freedom 00 . In an earlier paper uns, a Metric-Skewtensor-Tensor- 
Gravity (MSTG) theory was developed in which the effective coupling strengths G and y c of the gravitational 
field and a skew symmetric tensor field F were described by a a renormalization group flow formalism, 
analogous to the formalism in the standard model. A formulation of a gravity theory in which an effective 
average action T*. is a “course grained” Wilsonian action function with a variable momentum scale k has 
been developed by Reuter and collaborators Hung Ha ei ei use]. The theory was formulated to yield 
asymptotic safe quantum gravity removing divergencies in the ultraviolet (UV) |18| . 

In the following, we will adopt an effective running of the coupling constant G(fc), the gravitational 
charge Q(k) and the mass /i(k) in a version of the modified gravity (MOG) theory, assuming that there is 
an underlying quantum theory of the classical action. Such a microscopic interpretation of the theory is not 


1 



presently available, so we adopt a semi-phenomenological description of the scaling behavior of the MOG 
constants, which allows for an application of the model to the presently available early universe cosmology 
and late time astropliysical and astronomical data. This resolves the problem of the “missing” mass in the 
present universe, which is one of the most pressing problems in modern physics. As well as the metric g^, 
the repulsive gravitational vector field and its associated field strength B play a prominent role in the 
theory as dynamical fields. The gravitational coupling constants G and Q, the cosmological constant A and 
the vector field mass-range parameter g are described as scale dependent running constants, obtained by 
solving appropriate RG flow equations. 

In general relativity GR there is no known way to establish a Fourier transform linking momentum space 
to position space mmm A Fourier transform exists in fixed background symmetric spaces such as 
Minkowski spacetime and in Friedmann-Lemaitre-Robertson-Walker (FLRW) spacetime. The formal lack of 
a link between the particle physics momentum space and position space has been a source of conflict between 
GR and quantum gravity (QG) and this conflict can only be removed when a fully consistent QG is obtained. 
Making an analogy of the running of the gravitational coupling constants such as G(k ) with quantum field 
theory (QFT) and quantum chromodynamics (QCD) only applies in a semi-phenomenological framework. 
We describe the effective action T*, as a course-grained average description in the RG flow formalism. The 
effective course-grained approach advocated by Reuter and collaborators can provide a phenomenological 
RG flow treatment of gravity that can be expected to be valid in nature. 


2 The Gravitational Action and Field Equations 


The diffeomorphism invariant action for the theory is given by 

S = Sq + S(f, + S^m + Smi (1) 

where Sq is the Einstein-Hilbert action (we choose units such that c = 1): 

5g = \Lg J (]ix ^ {R + 2A )’ ( 2 ) 

where g = det(g Ml ,), g^ is the symmetric metric tensor of pseudo-Riemannian geometry, R = g^R^u, is the 
Ricci scalar and A is the cosmological constant. Moreover, Sq, is the field action: 

s <P = -^j > ( 3 ) 

where 

(4) 

We shall choose the potential V{<j>y) to be 

V{^) = -\n 2 r^ + W{^), (5) 

where W{<f>^) denotes a q field self-interaction potential. We can choose as a model of the self-interaction 
potential: 

w{^) = \h(r^) 2 , ( 6 ) 

where h is a coupling constant. 

The gravitational constant G in the action Sg is defined in terms of the “bare” Newtonian constant Gn- 

G = G n Z , (7) 


where Z = 1 + a corresponds to a “renormalization” of G and a is a dimensionless parameter. 
The actions 5^, S^m and Sm satisfy the relations 


1 SSm _ 1 m 1 _ 1 T SS^m 

7 =^^ -~2 Mmu ’ ~^gj^- 2 ^"’ ~S^~ 


( 8 ) 
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Here, Tm^u is the energy-momentum tensor for matter, while is the energy-momentum tensor containing 
the contributions from the massive vector field <f > M : 




1 :: 


B u a B v 


n dV (^) (1 B A 


Moreover, = (J°, J 1 ) (i = 1, 2, 3) is the current density for the field gravitational charge: 


J M = y/ uGnPmu 

where pm is the density of matter and = dx^/ds. We have 

Q = VuG n M = J d 3 xJ°. 

The field equations derived from the action principle are given by 

G [iv — SttGT^ , 

dV(^) 






= -J 11 -. 


and 


^irBfn, + V ^iB^u + — 0. 


(9) 

( 10 ) 

( 11 ) 

( 12 ) 

(13) 

(14) 


Here, G M „ = ^ g^R, T IJL , = Tm^ + and V a denotes the covariant derivative with respect to g 

The Bianchi identities satisfied by the Einstein tensor, \/ v G fJ ’ 1 ' = 0, lead to the conservation laws: 

V„T^ = 0. (15) 


3 Renormalization Group Flow and Running Coupling and Mass 
Constants 

From the point of view of QFT, we know that the electric charge e or the fine structure constant ctQED 
are not really constant but they are scale dependent or “running” quantities. In quantum electrodynamics 
(QED) the coupling constant olqed depends on the renormalization scale fc, a parameter with the dimension 
of momentum. The QED vacuum is a sea of virtual positron-electron pairs. A test charge is screened by 
a cloud of positron-electron pairs surrounding it, so that at smaller distances the test charge appears to be 
larger due to screening and smaller at large distances. In the non-Abelian Yang-Mills theory or quantum 
chromodynamics (QCD), the strong coupling constant a s = a s (k) is anti-screened, so that at large distances 
a s (k) grows with distance. In analogy with the situation in QCD, we can envisage that the gravitational 
coupling constant G = G{k ) is also anti-screened, because of the non-Abelian nature of the gravitational 
field. Then, G = G(k) is expected to grow in size as the distance scale increases due to non-perturbative 
renormalization effects. 

We have postulated a renormalized Newtonian gravitational constant G = G^Z in our action Eq.([2]): 

G = GZ = G N (l+a). (16) 

We shall now implement the running of the effective G and the B^ field gravitational coupling constant 
a Q = Q 2 /4tt [h = 1 and Q = \ZoGnM) by using RG flow arguments in which the “classical” coupling 
constants G and Q possess a scale dependent running behavior obtained by solving appropriate RG flow 
equations. 

We postulate scale dependent effective actions T^g^ and F* [0^], corresponding to “coarse-grained” 
free energy functionals, which define an effective field theory valid at the scale k or length scale £ = 1 /k, 
where £ is the length scale resolution of the observed system. Then, the Tare the bare actions obtained by 
integrating out all quantum fluctuations with momenta larger than an infrared (IR) cutoff kj r, or wavelengths 
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smaller than i. The T k are determined by functional differential equations for the RG flow, including the 
dimensionless coupling constants G(k) = k 2 G(k),aQ(k), the cosmological constant A (k) = A {k)/k 2 and the 
mass ft = n(k)/k. The dimensionless constant G(k)A(k) = G(k)X(k) is the quantity that is determined by 
measurement. We observe from CD and CD that both G(k) and ctQ(k) are determined by the running 
of a(k). The RG flow equations for a(fc), A(fc) and f{k) are described by the ordinary coupled differential 
equations: 

k-^a{k) = j3 a {a,\,fi), k—\(k)=/3x(a,\,p ,), k—f{k) = /3p,(a, A, ff). (17) 

To solve the problem, we are required to restrict the RG flow to a finite-dimensional subspace CD CD, 
corresponding to a truncated theory space. In the truncated space only the coupling constants G and ag, 
the cosmological constant A and the mass ft are considered. 

Every solution G(fc), ag(£:), A(fc), ft(k) of the truncated flow equations is associated with the family of 
action functionals: 

= Tfc^] +T k [(j) ll ]. (18) 

The effective average action formalism follows the RG flow from the bare action T^ =00 = S', corresponding 
to the initial condition for the Tfc-flow equation, down to k = 0 and T k~o = T corresponding to the effective 
action. We have to solve the effective equations of motion 

[dtiy] _ q STk[4>fj] _ Q Qg\ 

Sg^. v ’ Stpt* 

In practice, we truncate the space of solutions and for the UV sector, simple local truncations are sufficient 
to describe the physical system, but in the IR limit k —» 0 nonlocal terms must be included in the truncation 
procedure. 

We shall presently just consider the running of the coupling constants G{k ) and ag(fc). The RG flow 
equations are dominated by two fixed points G* and ag* : Gaussian fixed points at G* = ag* = 0 and non- 
Gaussian ones with G* > 0 and ag* > 0. The high-energy, short distance behavior is governed by the non- 
Gaussian fixed points, so that for k —> oo all the RG trajectories run into these fixed points by taking the UV 
cutoff along a trajectory running into the fixed points. This could lead to a non-perturbatively renormalizable 
quantum gravity theory. The conjectured existence of such a non-perturbatively renormalizable quantum 
gravity theory still has to be demonstrated. 

The running of the constants G(k) and ag is described by 

kdkG(k) = rjGG(a(k)), fcdfcag(fc) = 7ygag(a(fc)). (20) 

If tjg > 0 and r)Q > 0, then G(k) and oiQ(k) increase with increasing mass scale k (or decreasing length 
scale t = fc -1 as k -» oo), while they decrease for rjc < 0 and t\q < 0. In analogy with gauge theory, the 
first case refers to “screening”, while the second refers to “anti-screening”. Both the G(k) = Gn{ 1 + ce(k)) 
and ceQ(k) = a{k)G nM 2 /are determined by the one free parameter a(k). We adopt the ansatz that 
a(£ = k~ x ) has the behavior shown in Fig. 1. It is small at the scale of the solar system, so that G ~ Gn 
and Q ~ 0 which guarantees agreement with all the solar system experiments. It also guarantees agreement 
with the timing data for binary pulsars. 

As in QCD, we see that if we interpret t = 1/k as a distance scale, then both G(k ) and otQ(k) increase 
all the way from the UV energy scale to the IR energy scale as k —> 0. This corresponds to a quantum 
gravity anti-screening , playing the analogous role to the anti-screening in Yang-Mills QCD. This interpolates 
between the constant Newtonian value G(fci a b) ~ Gn at large k corresponding to short distances and the 
increasingly larger values of G(k) at larger distances for small k. We can as in QCD eliminate an explicit 
dependence of G(fc) and aQ(k) on a cutoff kc by choosing a renormalization or reference momentum mu. 
The equations for G(k) and ctQ^k) now contain only finite, physically measurable quantities depending on 
k and irtR. In our formulation of the RG flow, we do not necessarily assume that the asymptotic freedom 
condition G(fc) —> 0 as k —> oo and that a non-Gaussian fixed point exists. 
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Figure 1: A schematic depiction of the running of a{£) {£ = k 1 ). The dashed curve shows the behavior of 
a(£) for a{£) —> 0 as £ —> oo. 


4 Galaxy and Cluster Dynamics 


We expand g^ around the Minkowski spacetime metric g llv : 

g A 

where \h^, v \ « 1. The action for a test particle is 


S tp = - ( rri / dr + q / O,,^— ), 


dxJ‘ \ 
dr ) 


( 21 ) 


( 22 ) 


where m and q = \faG~Nm are the test particle mass and gravitational charge, respectively, and = (<j> o, <j>i). 
The spherically symmetric static solution for 4>o{r) is obtained from the equation = d(f>o/dr): 

K d— — A* 2 0o = 0. (23) 

r 

The solution is given by 

Mr) = -QSEfc/y (24) 

r 

where the charge Q = y/aG~NM. 

The equation of motion of a test particle is given by 


d 2 x> 1 dx a dx 13 

-d^ +T ^^ 


V dT ’ 


m 


(25) 


where r is the proper time along the particle trajectory. We assume that in the slow motion approximation 
dr/dr ~ dr/dt and 2 GM/r < 1, then for the radial acceleration of the test particle we get 


d 2 r GM qQ exp(—/xr) 

172 + - 2 “ =- 2 -v 1 + A tr ) 

at z m r z 


(26) 


For qQ/m = ciGnM and G = Gw(l + ct), the modified Newtonian acceleration law for a point particle is 
given by m®: 


a(r) = 


GnM 


1 Qt 


1 - exp(—r/r 0 ) 



(27) 
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where rg = l//x. The acceleration law can be extended to a distribution of matter: 


a(x) = — Gn J d' 


r3 ,p( x ')(x-x') 


r'|3 


[1 + a — aexp(—/u|x — x'|)(l -f /r|x — x'| 


( 28 ) 


The parameter a will run with distance scale in the RG flow scenario as shown in Fig. 1. At each distance 
scale k = l/£ the parameter a(k = l/£) is treated as a constant, so that at a given physical scale ^ p h ys the 
acceleration law (l27l) and (l28l) is valid. At the distance scale corresponding to the solar system, 4oiar ~ 40 AU, 
a(£) « 1, so that solar system experiments including the bending of light by the sun and the Cassini probe 
time delay observation agree with GR. Moreover, the binary pulsar timing observations, for distance scales 
corresponding to the compact pulsar dimensions and £< 8 kpc, will agree with the observations of binary 
pulsar timings for a (£) << 1. For galactic and cluster scales £q, fits to galaxy rotation curves and cluster 
dynamics are obtained with a = 8.89 ± 0.34 and p = 0.042 ± 0.004 kpc” 1 without dark matter GUSE]- 
Moreover, the bullet cluster data can be fitted without dark matter [F . 


5 Cosmology 


A description of early universe cosmology is obtained by adopting the homogeneous and isotropic FLRW 
background metric: 

ds 2 = dt 2 — a 2 (t) ^ ^ 1 ^ 2 + r 2 (d0 2 + sin 2 ddc/) 2 )^ , (29) 

where I\ = —1,0,1 (length” 2 ) for open, flat and closed universes, respectively. We use the energy-momentum 
tensor of a perfect fluid: 

Tfj, v = (p + p)u^,u v - pg (30) 

where is the total energy-momentum tensor in our RG flow framework, and p and p are the density and 
pressure, respectively. In the following, we assume a spatially flat universe K = 0. In the background FLRW 
spacetime, we have </>o ^ 0, </>j = 0 and = 0. Moreover, we have 

P = PM+P<j>, P = Pm+P <)>, (31) 


where pi and pi denote the components of density and pressure associated with the matter and the vector 
field. 

The Friedmann equations are given by 


2 8nGp : A 


H 2 = 


3’ 


a 47 tG A 

- =-^(p + 3p) + -, 

a 3 3 


where H = a I a. The energy conservation equation is 

din a 


p + 3-^-(p + p) = 0. 


(32) 

(33) 


(34) 


A linear perturbation on the FLRW background will connect the theory with anisotropy observations in 
the CMB. In the conformal metric with the time dp = dt/a(t ): 


ds 2 = a 2 (p)(dp 2 — dx 2 ), 


(35) 


the metric perturbations are defined in the Newtonian gauge: 

doo(x,t) = a 2 {t){l + 2$(x,f)), ffij(x,t) = a 2 (t)(l - (36) 


where >F denotes the gravitational potential. The perturbations of the vector field <f>^ are defined by 


<^(x,t) = a(t)(^ i (t) + <5c(v(x,t))(<M t ) = 0 ), 


(37) 
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In a purely baryonic early universe, the coupling of baryons to photons before and during the recombina¬ 
tion era will suffer damping, causing the collisional propagation of radiation from overdense to underdense 
regions. In order to obtain the expansion of the universe and structure growth to agree with observation, we 
invoke a cold dark matter in the form of the neutral field q , which plays the role of a ” dark photon” 7 ' with 
a sufficiently big mass The dark photon 7 ' interacts weakly with standard model particles through its 
gravitational strength coupling Q = \JaGNm^ and it has zero pressure. In our cosmology model we have 


fib 


8 TrG N p b 
3 H 2 


flej) 


87 T G N p$ _ 8 nG N p 

vac 

3# 2 ’ A = 


(38) 


where pb and p vac 
from neutrinos: 


denote the density of baryons and vacuum density, respectively. We also have a contribution 


_ 87 tGnPv 

~ 3 H 2 


(39) 


where denotes the density of neutrinos. We assume that the dark photon density dominates before 
the formation of the first galaxies, and the acceleration of the universe expansion is determined by the 
cosmological constant A. 

In Figure 2, we depict the running of a{a(t )) as the cosmic scale factor aft ) increases with time. We 
choose k = tf/aft) where £ is a constant and before the formation of the first galaxies, p</, > pb . We have 
for G = Gat(1 + a) that G ~ Gn and a < 1 before galaxy formation. After the epoch of galaxy formation 
Pb > P 4 > and the baryon density begins to dominate the density of matter. In Figure 3, we schematically 
depict the running of m^, with increasing aft), when it reaches the value = 2.6 x 10 -28 eV corresponding 
to the value p = 0.042 kpc -1 or 1/p = 24kpc used to fit galaxy rotation curves and the dynamics of galaxy 
clusters. 

The Newtonian potential in the early universe is given by 


= 47rGiv(pb(5b + pcdm<5cdm), 


(40) 


where k 2 denotes the square of the wave number, pcdm = P<t> and <5,; denotes the perturbation density con¬ 
trast for the baryon and CDM components. If p^ is sufficiently large, then = <5cdm will not be erased 
and the cosmological model yields agreement with data for the structure growth, rate of expansion of the 
universe, the CMB angular acoustical power spectrum. We show in Figure 4, the best fit MOG (A CDM 
model) TT power spectrum [ 22 ] . The formation of galaxies about 4-5 million years after the big bang occurs 
during the time when the modified gravity theory takes over with pb > p 4 , and galaxy dark matter haloes do 
not form. The data for the matter power spectrum determined by galaxy correlation functions can be fitted 
within modified gravity without dark matter [ 23 (. [ 2 ] . 


We have revived the original proposal for modified gravity that G, «q and p be considered as running 
parameters, as in RG flow JT0|. With this interpretation, the applicable values of the parameters in a 
particular situation depend on how the observed system is probed by telescopes and satellites, and the 
corresponding spatial resolution. Observations of a distant galaxy require a resolution of order kpcs, and 
the parameter values a and p will correspond to that scale. If one observes objects on Earth or in the solar 
system, then the spatial resolution will be of order AU or less. This is analogous to measurements of running 
coupling constants in particle physics, where the scale is the inverse of the momentum transfer, and higher 
energy collisions probe smaller scales. The geometry determined by the metric and curvature of spacetime 
is resolved spatially by observation and objects farther from observing instruments are not as well resolved 
as closer objects such as our solar system. 


6 Conclusions 

We have developed a gravity theory based on a modification of GR in which the action has in addition 
to the Einstein-Hilbert action for the metric g^ a dynamical action for a massive vector field and the 
field strength The gravitational coupling strengths G(k) and ctQ(k) run with momentum k in analogy 
with the running coupling constants in Yang-Mills theory and QCD. In our applications to astrophysics and 
cosmology, the G{k ) and ctQ(k) are anti-screened in the physical space of G, uq and A. The running of 
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increasing a 



of first galaxies 

a(t) 

Figure 2: Schematic running of a = a(a(t)), where a(t) denotes the cosmological scale factor. As a(t) —> 0 
a —> 0 according to the identification k = £/a(t) and G —tG^ in the Cosmic Microwave Background (CMB) 
and Big Bang Nucleosynthesis (BBN) epochs. 



a(t) 


Figure 3: Schematic running of the vector field mass mAa{t )) as a(t) —> oo. 
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Figure 4: The TT power spectrum using Q c h 2 = fl^h 2 , and the remaining parameters of the best 

fit listed in Table 4 of the Planck 2015 paper [22]. The red line shows the best-fit MOG spectrum (and the 
ACDM model spectrum). The lower panel shows the residuals with respect to the MOG model. 


the coupling constants G(fc) and otQ[k) is determined by Wilsonian RG trajectories. The trajectories have 
Gaussian fixed points when the fie and f3 aQ functions vanish at the points G* and oq* . Non-Gaussian fixed 
points can occur for non-zero values of G* and oq*. If they occur on the RG flow trajectories, it can be 
argued that unphysical divergences do not occur for UV values of k HU [HE]. The running of G(£) and of 
&q(£) is determined by the running of the dimensionless constant a(£), such that for small distance scales 
corresponding to the solar system or binary pulsars a < < 1. For cosmology the running of a is controlled by 
the magnitude of the cosmic scale aft) according to the scaling parameter k = £/aft ). In the early universe 
before the formation of the first galaxies a < 1, and the mass of the vector field is big enough to 
allow the vector particle be a hidden dark matter photon j', while for the late-time universe m$ ~ 10 -28 eV. 
In the late-time epoch the density p^ is diluted and p^ < p&, so that the baryon density p& dominates and 
the dark photon , coupled with gravitational strength to standard model particles, is undetectable. The 
rotation curves of galaxies and the cluster dynamics, including cluster collisions such as the Bullet Cluster, 
are fitted without dark matter mmm®. 

The RG flow gravity theory and the determination of the running of the parameter a are expected to 
be reliable for weak gravitational fields , while for the strong gravitational fields of black holes in the present 
universe mmm the running value of a will require a non-perturbative treatment outside the scope of 
the presently developed theory. 

For appropriate running of the constants G(fc), cxq and for physical values of the scaling parameter 
k or length scale £ = A; , observational data for the solar system, binary pulsars, galaxies, galaxy clusters 
and cosmology can be fitted by the theory. 
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